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We propose an efficient mechanism to heat a plasma by an intense microwave field solving the 
equation of ion motion in a wave field and a constant magnetic field in a large coupling regime. The 
mechanism does not relies explicitly on stochastic motion and is able to increase the ion velocity by 
several magnitude orders. Known thresholds for the onset of stochastic motion are also obtained. 

PACS numbers: 52.50.Sw, 52.20.Dq 

The possibility to realize confined nuclear fusion is linked to the ability to let ions reach enough high energies to 
ignite the process. Heating of such plasmas can be obtained by different means the most common of which is by 
radiofrequency directly injected through the confining system that is generally a tokamak. 

Such a heating exploits directly the interaction between a charged particle and an electromagnetic mode of a field 
propagating inside the magnetically confined plasma. This kind of interaction, being ruled by a strongly non-linear 
equation, is rather involved and a lot of analysis have been performed in the latest thirty years to properly understand 
the process. 

The most exploited effect that happens in such wave-particle interaction is the onset of a stochastic behavior 
according to Kolmogorov- Arnold-Moser (KAM) theorem [1 EH 0,0. 

The availability of a large portion of phase 
space to the ions, when a threshold is overcome, and the consequent entering into a higher energy regime is termed 
"stochastic heating" and it proves to be a very efficient method to heat a plasma by radiofrequency. 

Stochastic heating is commonly analyzed through a well-known model that, with the proper modifications, is able 
to give an adequate understanding of the physics at hand. Indeed, for an ion with charge q and mass m with a 
constant magnetic field B = (0, 0, Bq) and a wave E = (0, Eq, 0) cos(ky — cot), one has the equation^], 

y(t) + ^ 2 y(t) = —cos(ky(t)~Lut) (1) 

being Q = qB^/m® the cyclotron frequency. We make this equation adimensional by introducing y(t) — •> ky(t), t' = fit, 
a = Eok/BgCl and v = lj/Q so that, with the above redefinition of y, 

y(t')+y(t') = acos(y(t')-vt') (2) 

where we recognize that two critical parameters fix the physics for this problem, a and v. This is a Hamiltonian 
system. Indeed, it known that for stochastic heating produced by a lower hybrid wave the stochasticity threshold, 
assuming v > 1 , is given by 0, HI 
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while for ion cyclotron heating this condition becomes [| a « 1. Dissipation of the energy of the wave is due to 
Landau damping that, for the presence of the magnetic field, makes the absorption irreversible 0. 

The approach that is normally adopted to cope with this kind of problems is to work out canonical perturbation 
theory that holds for small a. Then, Chirikov overlap criterion for resonances and comparison with numerical 
results are used to determine the threshold of the onset of stochasticity to see when effective heating starts to set in. 
In this situation physics is fairly well-known and a diffusion equation can be derived describing the behavior of the 
plasma under such conditions. The overall mark of stochastic heating mechanism is the fact that the system under 
study is not integrable and KAM theorem does apply. 

Anyhow, we note that the opposite limit a — > oo gives us the interesting result that the system is ruled by the 
equation 

yo(?) = a cos(y (*')-<) (4) 
that is an integrable one. Indeed, the solution to this equation is straightforward to write down as 

7T 

Vo (t ) = 77 + vt' + 2 arcsin(msn(v / at', m)) (5) 
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being sn the Jacobi snoidal elliptic function and m its modulus. We notice that if we take 2/(0) = ir/2 and y(Q) = kvo/^l 
it easy to see that 
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(6) 
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and then the modulus of the Jacobi function depends on the initial velocity of the ion. Besides, T tr — 4_RT(m)/y / a, 
with K{m) = d<f>/y/T — m 2 sin (j), is the bounce time of a ion trapped inside a well of the wave. From this exact 
solution we can recognize two regimes. One has the ions moving uniformly with the phase velocity of the wave and 
the other is just the contribution of the trapped motion. If the first effect prevails one can have a net increase of the 
energy of the ion that, for the presence of the magnetic field, is irreversibly absorbed. Indeed, this can be assured by 
taking m <C 1 so that one can expand the snoidal function to have 

Vo(t') « ^ + id' + 2m sm(Vat') + 0(m 2 ) (7) 

and an efficient heating mechanism is granted. In fact, the averaged velocity on a bounce time will be given by 
{yo} = v - This means that the ion is moving with the phase velocity of the wave that is assumed much larger of the 
initial ion velocity. Similarly, one has {y 2 .) = v and a diffusion coefficient can be estimated dividing by the bounce 
time as D = {yl)/2T tr ~ ir^/y/a. 

This scenario would be consistent if we would be able to develop an asymptotic perturbation series in the limit 
a — ► oo to compute higher order corrections due to thepresence of the magnetic field. Indeed this can be accomplished 
by the duality principle in perturbation theory 0, [ij, [ToL [ill IT^ . This principle can be stated by saying that 
interchanging the perturbation terms in an equation gives two perturbation series having a development parameter 
that in a series is the inverse of the other. That is, in our case one would have (weak perturbation, a — > 0) 
y = yo + ayi + a 2 y-2 + ... and for the dual series y = yo + ot y% + a~ 2 y 2 + .. . that holds in the opposite limit a — •> oo. 
In order to get a non trivial set of perturbation equations, one need to rescale time because, as already seen 

above, this is the proper scaling of the leading order solution. So, we take 

t = VSf (8) 
Vo = yo + ot~ 1 y 1 + a~ 2 y 2 + ■ ■ ■ 

into the motion equation and the following working set of equations is easily obtained 

yo{r) -cos [ yo(r) -j=T ) = (9) 
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j/i (t) + sin ( t/o (t) - -j=t J yi (r) = -y (r) 

'") - sin ( 2/o(t) - -4=rj J/2 (t) = -t/i(r) - -cos fj/o(r) - -^=r ) //[(V) 



Using the exact solution to the leading order (J5J the above set of equations takes the form 

27i(r) + [l-2m 2 sn 2 (r,m)]2/i = -yo(j) (10) 
2/2 (t) + [1 — 2m 2 sn 2 (r, m)]j/2 = — 2/i( r ) + wsn(r, m)yl — m 2 sn 2 (T, m)y\{T) 

that in the limit m <C 1 take the simpler form 

yi{r)+yi = -Vo{t) (11) 
J/2(t)+2/2 = -Vi(r) +msin(r)2/ 2 (r) 



to the first order in m. These equations can be solved straightforwardly taking at the leading order eq.J7|). 
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The solutions contain secular terms. These terms can be removed in many ways. The one we prefer is the 
renormalization group method as presented in [T3L IbH Il5| . This technique assumes that the perturbation solution 
should be computed to an initial time t as to have y(r, r ) = yo( T ! r o) + a ~ 1 Ui( T J T o) + a2 U2( T , T o) + • ■ ■ and the 
initial constants depend on the initial time, <fi = 4>(r ) and m = m(r ). The requirement that the perturbation series 
should not depend on the choice of the initial point translates into the condition 



dy(r, t ) 



dr 



= (12) 



that gives back renormalization group equations for the constants <p and m that now turn out to depend on r. Then, 
the condition y(r, To)| TO _ T produces the correct perturbation series without secular terms. It is interesting to point 
out that all we have done is just to compute an envelope |l5j. 
After resummation we get 

Vo(t) = J + ^=r + 2m(T)sm(r + ^(r)) (13) 

Vi( T ) = |[cos(r + ^(r)) - 1] -m(r)sin(r + ^(r))) 

7T ( -K 2 3\ 

2/2(t) = -[1-cos(t + 0(t))]+ f — + -J m(r)sin(T + 0(r)) 

7T 2 7T 2 

— — m(r) sm(r + </>(t)) cos 2 (t + 4>(t)) + — ra(r) sin(T + <^>(t)) cos(r + </>(r)) 

where now the modulus m depends on r and the phase 4>{t) contains a frequency shift as from the renormalization 
group equations 

dm 7r „ / ^ N , , 



d4> _ l 1 /5tH 
dr ~~ 2a 4a 2 
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Having the modulus depending linearly on t implies that if one waits a time long enough our initial approximation 
m <C 1 in absolute value may fail. We discuss this point below about stochasticity thresholds. We have also found 
higher order corrections to the bounce time. Finally, we note the renormalization constant 1 — 1/a + I /a 2 + . . . — 
a/(l + a) to the it/2 term. 

We now discuss the stochasticity threshold in the case of a ion cyclotron wave (ICW) [l| and lower hybrid wave 
(LHW)0, Q- The mechanism of heating does work only if the particle moves freely at the phase velocity of the wave. 
Instead, if bouncing of the ion in a well of the wave prevails, heating is not effective. For a ICW one has v = 1 and 
the only guaranty of a fully consistency of the above discussion is a 3> 1 in agreement with the result given in 0. 
For a LHW, when a is not so large but greater than one and v ^> 1 as discussed in 0,0, it is fundamental that the 
modulus m does not prevail on the free particle term forcing the ion behavior to a simple bouncing in the wells of the 
wave. This is granted taking 1/v » Tr/8a^ giving us the improved threshold for this wave 

(7T2/)i 

a » { -^- (15) 

and the numerical factor is now about 0.54. In this way we have proved fully consistency of our approach with respect 
to previous ones on stochastic heating obtaining a deeper understanding of the involved physics. 

In order to complete our analysis we give here some numerical studies that are essential to support the above 
scenario. We have considered two cases for a, v and y(0). We have not taken a too large because otherwise the plots 
would not be much interesting as the analytical solution hits the exact one as does the numerical one. The goodness 
of the approximation we have applied is evident provided that m is kept small and the thresholds given above are 
sufficiently overcome. So, figH corresponds to have \m\ « 0.45 but a — 10 is not so large. The same happens in fig[21 
that has \m\ « 0.022. The worst situation is seen in figl^where |m| 1=3 2.3 and the thresholds not properly overcome. 
So due to the value of a — 10 and the corresponding problems in m and thresholds the agreement goes from relatively 
good to bad. The relevance of the threshold is due both to the largeness of a and to avoid that m increasing in 
time may overcome the initial condition m« 1 from which we started. Things improve, as should be, with a = 100 
where the agreement is very satisfactory. An interesting situation is seen in figEl where the solution given in eq.(|5|l 
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FIG. 1: Numerical, analytical and eq.lJ^J for a = 10, v = 15 and y(0) = 0.01. 
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FIG. 2: Numerical, analytical and eq.Q for a = 10, v = 0.15 and y(0) = 0.01.. 



is better than the perturbation series. This is due to the fact that, although the thresholds are properly overcome, 
one has m ss 0.75 very near one, where our fundamental assumption to derive the perturbation series starts to fail. 
As expected, in some cases the perturbation series performs much better than eq.©. 

So, we can conclude that a satisfactory perturbation approach has been presented here to treat problem of strong 
particle- wave interaction in a plasma. A heating scenario for ions has been given showing that it can increase the ion 
velocity of several magnitude orders making it move to the phase velocity of the wave. Known thresholds have been 
derived in a fully analytical and rigorous way forming a complete scenario for plasma heating by radiofrequency. A 
lot of problems having a longstanding numerical treatment could be faced now in an analytical way by the method of 
dual perturbation theory giving an improved insight into physics. 
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FIG. 3: Numerical, analytical and eq.0 for a = 10, v = 0.15 and y(0) = 3. 
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FIG. 4: Numerical, analytical and eq.@ for a = 100, v = 15 and y(0) = 0.01. 
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